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ABSTRACT 

Ram pressure stripping is an important process in the evolution of both dwarf galaxies and large 
spirals. Large spirals are severely stripped in rich clusters and may be mildly stripped in groups. 
Dwarf galaxies can be severely stripped in both clusters and groups. A model is developed that 
describes the stripping of a satellite galaxy's outer H I disk and hot galactic halo. The model can be 
applied to a wide range of environments and satellite galaxy masses. Whether ram pressure stripping 
of the outer disk or hot galactic halo occurs is found to depend primarily on the ratio of the satellite 
galaxy mass to the mass of the host group or cluster. How the effectiveness of ram pressure stripping 
depends on the density of the inter-group gas, the dark matter halo concentrations, and the scale 
lengths and masses of the satellite components is explored. The predictions of the model are shown 
to be well matched to H I observations of spirals in a sample of nearby clusters. The model is used to 
predict the range of H I gas fractions a satellite of mass M v<sa t can lose orbiting in a cluster of mass 
M V g r . 

Subject headings: galaxies: evolution — galaxies:clusters:general — galaxies: dwarf 



1. INTRODUCTION 

Ram pressure s tripping was first proposed by 
iGunn fc Gotti l)1972f ) to explain the observed absence of 
gas-rich galaxies in clusters. They noted that galaxies 
falling into clusters feel an intracluster medium (ICM) 
wind. If this wind can overcome the gravitational at- 
traction between the stellar and gas disks, then the gas 
disk will be blown away. They introduced a simple ana- 
lytical condition to determine when gas is lost: 

PiCMV 2 sat < 2TrGa*<j gas . (1) 

The right-hand side, 2irGa*<j gas , is a gravitational restor- 
ing pressure, picM is the density of the ICM, and v sa t is 
the orbital speed of the satellite galaxy. Using this con- 
dition they concluded that spirals should lose their gas 
disks when they pass through the centers of clusters. 

Galaxies near cluster centers have earlier type mor- 
phologies, are redder, and form few er stars than galax- 
ies in the field of s imilar luminosity iGpm gz et al1l2003t 
IGoto et a.U 12001: IHogg et a.U l20M 120(141) 1 Stripped 
galaxies have low H I masses for their size and mor- 
phology. Such H I defici ency correlates with lower star 
formation rates (SFRs) ijGavazzi et al.l l2006L and low- 
ered SFRs lead to redder colors. Therefore, ram pressure 
stripping, in conjunction with processes like harassment 
that can alter the structure of a galaxy, may play a role 
in transforming galaxy morphologies. 

Several sets of observatio ns indicate that stripp i ng oc- 
curs to spirals in clusters. iGiovanelli fe Havnesl l|1983f) 
define a deficiency parameter that compares a galaxy's 
observed H I mass to the expected H I mass in a 
field galaxy with the sa me morphology and optica l 
size. Galaxies in Virgo. l|Giovanelli fe Havnes) [1983L 
Com a. IIB ravo-Alfaro et al]l2000j) , and other nearby clus- 
ters. llSo lanes et alJ2001l) are observed to be H i deficient. 
iSolanes et al.l l)2001|) observe that the average deficiency 
increases with decreasing distance to the clusters' centers 
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and that de ficient galaxies have, on average, more eccen- 
tric orbits. iCavatte et al.l p994) observe indications of 
stripping in the H I distributions of individual galaxies 
in the Virgo cluster. Compared to an average profile for 
a field galaxy of similar morphology and optical extent, 
these galaxies have normal gas densities in their inner 
disks but are sharply truncated. More detailed obser- 
vations of several Virgo spirals have since been carried 
out. NGC 4522 is observed to have an undisturbed opti- 
cal disk, a severely truncated H I disk, and extra planar 
Ha and H I on one side of t he disk ijKennev fc Koopmanl 
Ii999t iKennev et all 12004^ . IC 3392, NGC 4402, NGC 
4419, and NGC 4388 all have a trunca ted H I disk 
and th e first three have extraplanar H I iJKennev et alJ 
l2004a|K NG C 4569 has a trunca ted H I disk and extra 
planar Ha ijKennev et alJl2004ajK Observing asymmet- 
ric, extraplanar gas in combination with an undisturbed 
optical disk is a strong indication that the galaxy's gas 
is interacting with the ICM. 

Detailed simulations of ram pressure stripping con- 
firm that spirals in clusters should be stripped. For 
some simulations the final radii of the galaxies are com- 
pared to the predictions of the Gunn and Gott condi- 
tion. In addition, the morphologies of simulated galax- 
ies have been compare d to observa t ions o f spirals in 
Virgo. Simulations by lAbadi et al.l l)1999j) show that 
galaxies in a face-on wind are stripped to radii near 
those predicted by the Gunn and Gott condition, but 
galaxies in an edge-on wind are only mildly stripped. 
In contrast, other simulations show a two-step strip- 
ping process that strips galaxies at all inclination an- 
gles. Ram pressure stripping first quickly strips gas from 
the outer disk, and then the remaining d i sk is slowly 
viscously stripped iMarcolini et al.l l2003t iQuilis et alJ 
200(1 IRoediger fc Henslerl 120051 ISchulz fc Struckl l200D . 
Marcolini et al.l ((2003) point out that, unlike ram pres- 
sure stripping, viscous stripping is more effective in an 
edge-on wind. They find that the final radii of the face- 
on galaxies match those predicted using the Gunn and 
Gott condition and that edge-on galaxies are stripped 
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to slightly larger radii. iRoediger fc Henslerl ((2005) run a 
series of simulations to test the dependence of the strip- 
ping radius on the speed of the ICM wind, the mach 
number, and the vertical structure of the gas disk. They 
find that the most important parameters are the ram 
pressure and surface density of the gas disk and that the 
effects of varying the mach number and vertical struc- 
ture are minor. They find that the original Gunn and 
Gott condition does a fair job predicting the stripping 
radius. However, when they use an adjusted ram pres- 
sure, they find that using the thermodynamic pressure 
in the central plane of the disk more accurately predicts 
the stripping radius. 

Simulations of ram pressure stripping can also 
match the range of morph ologies of stripped galaxies. 
IVollmer et"aH ij2000L I2004J) track the density, velocity, 
and velocity dispersion of a simulated galaxy's cold disk 
gas. In NGC 4522, the kinetically continuous extrapla- 
nar gas matches simulations of a galaxy that is currently 
being stripped. However, the gas in NGC 4569 is less 
pronounced, not kinetically continuous with the disk gas, 
and has a large velocity dispersion, more like simula- 
tions of a galaxy that has been st r ipped in the past. In 
their simulations iSchulz fc Struckl 1)200 lj) observe a pro- 
cess they term "annealing". The ICM wind compresses 
the gas disk and triggers the formation of large spiral 
arms. The interaction of these arms with the wind leads 
to the contraction of the galaxy and the formation of a 
dense truncated gas disk. A burst of star formation may 
occur at the edge of the truncated disk. A ring of star 
formation is seen at the truncation radius of the disk in 
IC 3392. 

Some groups have studied stripping in groups and clus- 
ter outskirts where the ram pressure is lower than in 
cluster centers. While simulations and observations both 
support the idea that large spirals in clusters are stripped 
of their H I disks, it is not clear that stripping occurs 
in poorer environments. Two types of stripping may 
occur in groups. The ram pressure in groups is lower 
than in clusters. However, dwarf galaxies have lower 
masses than large spirals and therefore lower restoring 
pressures, and dwarf galaxies may be stripped of most 
of their H i. In addition, large spirals have restoring 
pressures that decre ase with radius and may be stripped 
of their outer disks. iMarcolini et al ] ll200l run simula- 
tions of dwarf-like galaxies in winds typical of the ram 
pressure in groups and find that the simulate d dwarfs 
are stripped. The IRoediger fc Henslerl l|2005D simula- 
tions include runs with low ram pressure. Their large 
and medium-sized model spirals are stripped of some of 
their outer H 1 in these runs. The observational evidence 
for stripping in low-mass systems is circumstantial. The 
gas content and star formation histories of Local Group 
dwarfs correlate with their distances from th e dominant 
spirals IBlitz fc Rohisha,wll27MlGTeT^l27iril. 

This project aims to use the clear observations of 
stripping in clusters to predict when stripping occurs in 
groups. This is done by examining how the ram and 
restoring pressures depend on the masses of the cluster 
or group and the satellite galaxy. Two ways in which the 
masses of the host system and satellite galaxy affect the 
ram and restoring pressures are identified. Satellites' or- 
bital speeds are set by the depth of the cluster's potential, 
and the restoring force on the gas disk is determined by 



the depth of the satellite's potential. In addition, both 
cluster and galaxy morphologies change with mass. A 
model is developed and used to study both types of mass 
dependence. This is done in the following steps. 

In the next section a model for the mass and gas distri- 
butions is presented both in terms of physical parameters 
and in terms of scale-free parameters. In section 3 the 
dependence of ram pressure stripping on the masses of 
the satellite and cluster is presented assuming that the 
scale-free parameters are constant. In section 4 the de- 
pendence of the scale-free parameters on the masses of 
the satellite and cluster and the sensitivity of the model 
to these changes are explored. In section 5 the results 
are discussed and the model is compared to observations 
of stripping. Section 6 concludes. 

2. THE MODEL 

The analytical model of ram pressure stripping devel- 
oped here uses the Gunn and Gott condition to predict 
the extent to which galaxies are stripped. This condition 
is based on a simplified picture of stripping that ignores 
the details of the hydrodynamical processes and balances 
forces on large scales. However, simulations of galax- 
ies that are stripped by a face-on wind result in strip- 
ping radii that match the predictions of the Gunn and 
Gott condition llAbadi et alJll99ft IMarcolini et~aTll2003t 
IRoediger fc Henslerl 120051) In the simulations galaxies 
in an edge on wind are eventually viscously stripped to 
radii similar to those of face-on, ram pressure stripped 
galaxies. This is discussed further in section EH I n light 
of the general agreement between Gunn and Gott's con- 
dition and simulations of ram pressure stripping, it is as- 
sumed that a model based on this condition can be used 
to make general statements about where and to what 
extent galaxies are stripped. 

The model for the group or cluster has two compo- 
nents, a gravitational potential in which the satellite or- 
bits and an ICM against which the satellite shocks. The 
gravitational potential is modeled using an NFW profile 
with mass M Vtgr , radius, r„ i9r , and concentration, c g . 
The scaled radius s = r/r Vigr . 

-GM VtBT ln(l + c g s) 

0halo = —g{Cg) — (2) 

Tv,gr ° 

The NFW profile is discussed in more detail in ap- 
pendix E| The ICM is modeled using a /3 profile with 
scale length r c , and central density, pq. 

/ 3 \ (-3/2)/3 

P = Po ( 1 + ^ J ( 3 ) 
The (3 profile can be rewritten as 

p =2Mc( 1+ 3 (4) 
3 V {re/r v ,g r r ) 

The parameter a relates po to the characteristic density 
of the NFW profile, and v p p® is the average dark matter 
density within r Vtgr . 

The gravitational potential of the satellite and the sur- 
face density of the gas disk are needed to model the 
restoring pressure. The gravitational potential combines 
a dark matter halo, a stellar disk, and a stellar bulge. 
The dark matter halo is modeled as an NFW profile 
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with mass M VtSat and concentration c s . The disk, of 
mass Mrf, is described using a M iyamoto and Nagai po- 
tential {Miyamoto fc Nagailll975|) . 



-GM d 



(5) 



R? + (R d + Vz 2 + h ,2 Y 



This potential is easily differentiated, and with the ap- 
propriate choice of disk scale length, R d , and height, h, 
many of the properties of the potential of an ex ponen- 
tial disk can be matched {Johnston et alJll9 95). This 
potential can be rewritten as 



-GM d 



1 



™disk 



(6) 



where X d = r V:Sat /R d , X h = h/R d , S = R/r VjSat , and 

%s — Z /fv,sat- 



The bul ge, of mass M b , is described using a Hernquist 

-GM b 



a,s s ivih,y t 

potential {Hernauistl ll990) . 

(fibulge - 



r + n 



This can be written as 



-GMh 



1 



n hs S at + 1 



(7) 



(8) 



where X b = r VySat /r b , s sat = r sat /r Vl and r sat is the dis- 
tance from the center of the satellite. By introducing 
m ds = Md/M ViSat and m b = M b /M VySat , the full grav- 
itational potential of the satellite can be written in the 
form 

i GM V sa t £ / A \ /r .\ 

9 full = J [Ssat, mi, Xi,c s ) (9) 

Tv ,sat 

Gas in galaxies can be found in three spatial compo- 
nents: an exponential molecular gas disk with a scale 
length comparable to the stellar disk, a nearly flat atomic 
disk that extends beyond the stellar disk and shows a 
sharp cutoff, and a hot galactic halo. Beyond the disk 
cutoff, the atomic disk may co ntinue as an ionized gas 
disk {Binnev & Merrifieldlll99^ . 

The H I is modeled as a thin flat disk with mass M g , 
surface density, a g , and a sharp cutoff at radius R g 



M g 
TrJZf 



(T , R < Rg 
0, R > Rg 

_ M v , sat m dg X 2 g 



(10) 



(11) 



where m dg = M g /M v sat is the fractional mass of the gas 
and A g = r v<sa t/ R g is the scaled size of the disk. 

The inner gas disk is dominated by the H2 disk. The 
molecular disk is significantly more difficult to strip than 
the H 1 disk both because the H2 disk is more compact 
and because the H2 is found in molecular clouds. These 
clouds do not feel the effect of stripping as strongly as 
the diffuse H 1. How stripping of H 2 occurs and how 
the H2 clouds and diffuse H 1 in the inner disk interact 
are not known. The two phases may be tied together 
by magnetic fields; in which case the inner H 1 disk will 
remain until the ram pressure can remove the entire inner 
disk. It may also be possible for the wind to remove the 



H 1 from the inner disk while leaving the H2 behind. 
Because of this uncertainty, the stripping of the inner 
disk will not be studied here. The model is only used to 
discuss stripping beyond 1.5 stellar scale lengths. If the 
ram pressure is not strong enough to affect the H2, the 
H2 will contribute to the restoring pressure in the outer 
disk. Therefore, the mass of the H2 disk is added to the 
stellar mass. 

The emphasis of this project is to study the stripping 
of the H 1 disk. However, the fate of the hot galactic halo 
can also influence the evolution of the stripped galaxy. It 
is expected that the hot halo will be easily stripped both 
because it is diffuse and because the restoring potential 
of the satellite is strongest in the disk. To check this a 
h ot halo is included in th e model galaxy. 

iMori fc Burkertl l)2000D give an analytic condition for 
the complete stripping of a hot galactic halo. 



2 ^ u P0,satkBT 
PlCMV sat > f o t h = 



(12) 



where Po,th and po,sat are the central thermal pressure 
and density of the hot galactic halo, ks is the Boltzmann 
constant, T is the temperature of the galactic halo, and 
\xm v is the average molecular mass. The thermal pres- 
sure replaces the gravitational restoring pressure in the 
original stripping condition. In their simulations, this 
condition does a fair job of predicting the mass of a 
galaxy that can be completely stripped. This condition 
is adapted for the current model and all gas outside the 
radius where P ra m = Pth is assumed to be stripped. 

The hot galactic halo is modeled by assuming that the 
gas is at the virial temperature, T v , of the satellite and is 
in hydrostatic equilibrium in the satellite's dark matter 
potential. The self-gravity of the gas is ignored. The 
density of the hot halo is 



P(Ssat) 



0P0 



m S gj(s sa t, Cs) 



(13) 



where scaled density profile j(s sa t,c s ) is defined in ap- 
pendix [C] and m sg is the ratio of the mass of gas in the 
hot galactic halo to M VySat . 

A dimcnsionless temperature, t(c s ) is also introduced 
in appendix [HI It is defined by 



t(c s ) 



p/rrip 

It is a function of only c s . 



\GM V , S 



(14) 



3. RESULTS 

In this section the dependence of stripping on the 
masses of the large and satellite halos and on the model's 
parameters is explicitly determined. 

The maximum ram pressure, the ram pressure at the 
pericenter of the galaxy's orbit, is given by 

2 e av pP » ' ~ 2 N ( - 3/2)/3 



P, 



ram,max 



v ,9 r 



1 



( 15 ) 

where so is the distance of closest approach, e is the or- 
bital energy per unit mass, e v = -GM, ijr /r, !9r , b s = 
l 2 /(\e\ry gr ), and I is the orbital angular momentum per 
unit mass. Both b s and e/e v are expected to be indepen- 
dent of the masses of the satellite and the cluster. 
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At any point in the orbit, s or bit, the ram pressure is 
2GM v , gr av p p° c 



-Pram y^orbit) 



P\Sorbit i Cg 



(16) 



where dimensionless pressure, p(s,c g ), is defined. 



p(s,c g 



ln(l + c g s) 



(r c /r v 



Equations for the orbital speeds are derived in ap- 
pendix El Assuming that a, 0, r c /r Vygr , and c g do not 
depend on the mass of the halo, the M v>gr dependence 

Of both Pram (So) and P ram {s orbit) is Pram OC Mv'gr- 

For a potential of the form in cq. [§1 the force per unit 
gas mass in the z direction is found as follows. 

GM V sa t d . _ 
Iz = ttJ (s sa t,rni,Xi,c s ) 



Assuming that this set of parameters is constant is 
equivalent to assuming that the physical parameters scale 
with mass in the most obvious way. It assumes that 
component masses scale as M v , lengths scale and 
the central densities of the ICM and the satellite's hot 
galactic halo are constant. For any model in which these 
(-3/2)/? scalings hold, the result that the extent of stripping de- 
pends on M VySat /M Vygr will hold. For any cluster mass, 
the scaled radius will set the ICM density and the or- 
bital velocity will be proportional to M Vigr /r Vigr . The 
restoring pressure depends on the depth of the satellite 
galaxy's gravitational potential well and the density of 
the gas disk. For a generic potential, oc M v ^ at /r v , sat 
and the restoring force per unit gas mass is proportional 



(17) 



to M. 



fz 



r v , sa t dz 
GM VtSat d 



(18) 



(19) 



Assuming that the m.j and A^ introduced in section|21and 
c s are constant with mass, the mass dependence of the 

1 /3 

restoring force per unit gas mass is given by f z oc M v ' sat . 

For the gas disk, if the fractional mass of the gas, m dg , 
and the scaled size of the disk, A g , are both constant, 
then 

t/r 2 cx Ml'l t (20) 



Og (X M V>Sat / 



Combining the above, 

Prest,max — &g fz,max OC ^v,sat (^-0 

For any radius along the disk, R s t r , there is a maximum 
restoring pressure. If the maximum restoring pressure is 
greater than the ram pressure, the satellite holds the gas 
at this radius. The condition for the satellite holding its 
gas can be written as 



rest,max 



M 2/3 (si") f {s S at,m ll \ i ,c s )2rn dg \ 2 g 



p, 



ram.max 



M 2/3 



> 



Otp(s or bit, C g ) 

ap{s orbi t ,c g ) 



> 1 



(dr) f {ssat-,rn il \i,c s )2m dg \ 2 g 



(22) 
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The restoring pressure of the galactic hot halo is 



(23) 



Pre. ftt — 



t(c s )m eg j(s sat ,c s ) (24) 



Tv,sat J \ 3 

The condition for retaining this gas within s sa t is 



M, 



v,sat 



> 



v.gr 



t(c s )m sg j(s sat ,c s ) 

2ap(s orbit, c g ) 
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(25) 



If the model's parameters: c g , c s , e/e v , b s , (3, r c /r Vigr , 
a, the rrii, and the A^; are all independent of the mass 
of the satellite and cluster, then the fraction of gas 
that is stripped from a satellite depends on the ratio 
Mv,sat/Mv,gr- In this case, P ram oc My( gr and for both 

3/2 

the gas disk and hot galactic halo P re st oc M v sat . 
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If the mass of the gas disk scales with M v , 



and the disk length with r v<sa t, then the density of the 
gas disk scales as M\^ at and the restoring pressure as 



M. 



2/3 



4. MODEL PARAMETERS 



In this section a reference set of values for the model's 
parameters is introduced and the effect of varying these 
parameters is discussed. The reference parameters de- 
termine the values of M VySat / M v>gr at which stripping 
occurs. They are taken from observations, where pos- 
sible, and ACDM simulations, otherwise. Two types 
of parameter variations are discussed. First, several of 
the model parameters vary systematically with M v>sa t 
or M Vj g r . Second, at fixed masses the parameters have 
scatter. 

The scatter in the parameters is discussed for several 
reasons. By combining an estimate for the scatter in each 
parameter with the model, the parameters that have the 
most effect on the radius to which an individual galaxy 
is stripped can be identified. Second, it is useful to know 
the range of ratios across which stripping occurs. Finally, 
if the overall scatter is not too large, the relationship be- 
tween R str and M VtSa t/M Vigr could be used to search for 
signs of stripping in a large galaxy survey. The results 
of this section are discussed in the terms above in sec- 
tion 

In Figures ^ an d the M V:Sat at which a satellite is 
stripped is plotted versus R s t r /R g and M str /M g for a 
satellite orbiting at a variety of s or bit values in 10 15 , 
10 14 , and 10 13 M Q clusters. The mass of the stripped 
gas is denoted as M str - Table 1 lists the model pa- 
rameters used for three models, a large satellite in a 
10 15 M© cluster, the reference model; a large satellite in 
a 10 14 Mq cluster, the middle-mass model; and a small 
satellite in a 10 13 Af Q group, the low-mass model. The 
large satellite model is based on the Milky Way. Physical 
lengt h scales and masses for the stellar disk and bulge are 
from lJohnston et aD (^95). The stellar disk is modeled 
using M ds = 10 n M Q , R d = 6.4 kpc, and h = .26 kpc. 
The equivalent R d for an exponential disk is 5 kpc. The 
bulge is modeled using Mb = 3.4 x 10 10 M Q and r\, = .7 
kpc. A flat disk with a radius of 30 kpc gives R g /R d 
similar to that seen in observations. The H I gas mass 
is chosen to result in a disk density of 5 M@ pc~ 2 . The 
virial mass is chosen such that M VjSa t/Mb aryon ~ 12. 
Setting M ViSa t also sets r VySat . The concentrations used 
are c s = 12 and c „ = 5.5, 7.5, a n d 10. The concentra- 
tions are based on iBullock et alj <)2001|) and assume an 
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Fig. 1. — Mass at which a satellite's gas disk is ram pressure stripped versus Rstr/Rd- Short-dashed line: s or bi t = .25, solid: s or tit = -35, 
long-dash: s or tit = .5, dash-dot: s or tit = .75, dash-triple-dot: s or tit = 1- Top-left: Spiral galaxy orbiting in a 1O 15 M0 cluster. Top-right: 
Spiral galaxy orbiting in a 10 14 Mq cluster. Bottom-left: Spiral galaxy orbiting in a 10 13 Mq group. Bottom-right: Dwarf galaxy orbiting 
in a 10 13 Afo group. 



average over-density within r v , v p , of 337. The assumed 
over-density affects the concentrations, which are mea- 
sured using v p — 337, and a, because it relates the ICM 
and dark matter densities. Note that v p does not appear 
in equation 1231 or f?ol The ICM parameters, a, (3, and 
r v,gr/ r ci are set by X-ray observations of clusters and 
groups, and the orbital parameters are taken from simu- 
lations. Both are discussed below. Table 1 also lists the 
assumed scatter in each model parameter and resulting 
scatter in the mass at which a satellite is stripped. 

4.1. ICM Profile 

The profile of the ICM is determined by the parameters 
a, r v ^ gr /r c , and (3 (eq. 0. Sanderson fc Ponmanl (|2003f) 
compile a set of average ICM profiles for clusters with 
temperatures ranging from 0.3 to 17 keV. They find that 
the ICM of groups and poor clusters is bot h less dense 
and m ore extended than in rich clusters. IMohr et alJ 
(1999) present X-ray observat ions of a set of clusters 
and | Helsdon fc Ponmanl l|2000j) and lOsmond fc Ponmanl 
l)2004|) present X-rav observations of groups. The three 
groups present (3 and r c (in units of kpc) for a best- 
fit j3 profile and the X-ray temperature, Tx, for each 
group or cluster. The average (3 and r v _ gr /r c for the 
cluster observations are approximately 0.67 and 12. The 
averages for the group observations are approximately 
0.45 and 170. To convert r c (in kpc) to r c /r Vigr , virial 
masses for the groups and c lusters are found usi ng the 
M 2 oo (Tx) relationship from iPopesso et all ((2005). The 
mass Af2oo is the mass contained within the radius at 
which the average over-density reaches 200. It is con- 



ve rted to M„ using equa tions IA7I and IA6I The profiles 
in lSanderson fc Po nman| (|2QQ^Ij are sc aled using r2oo as 
determined in ISanderson et alJ i|2003l) . The differences 
in the r c /r v , gr determined using the two Mzoo(Tx) re- 
lationships are within the intergroup scatter in r c /r Vtgr . 
The average profiles for the hotte st and coolest groups 
from ISanderson fc Ponmanl l)2003(l are well matched by 
using a = 6.5 and 20 respectively and the appropriate 
(3 and r v>gr /r c . In Figures ^ and El plots are made for 
stripping in 10 13 , 10 14 , and 10 15 M Q clusters. A cluster 
with M„ nr = lO 15 -/^ falls into the h ottest temperature 
bin from lSanderson fc Ponmanl l|2003|) while a group with 
M v>gr = 1O 13 M falls into the coolest. The ICM param- 
eters for the 10 14 Mq cluster are a — 5.5, (3 — .57, and 

r v,grl^c — 25. 

Scatter in these parameters is estimated using the three 
sets of observations. The observed scatter in (3 is ±0.08. 
The majority of clusters have r Vtgr /r c between 7 and 20, 
and groups have r v ^ gr /r c between 75 and 270. For both 
groups and clusters, a is assumed to vary by a factor of 
V6 in either direction. Varying (3 varies the ram pressure 
at Sorbit — .35 by a factor of 2 for the clusters and a factor 
of 7.4 for the groups. The ram pressure varies by a factor 
of 6 in the clusters and 5.4 in the groups when r v gr jr c is 
varied. Varying a by a factor of 6 results in an increase 
in the ram pressure by the same factor. 

4.2. Orbits 

The orbital parameters are important for determining 
both satellites' pericenters and orbital speeds. In theory, 
the pericenter of a satellite's orbit is determined by e/e v 
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Fig. 2. — Mass at which a satellite's disk is ram pressure stripped versus the mass fraction of gas disk that is stripped. Lines are the 
same as in Fig. 1. Top-left: Spiral galaxy orbiting in a 10 15 Mq cluster. Top right: Spiral galaxy orbiting in a 10 14 M o cluster. Bottom 
left: Spiral galaxy orbiting in a 10 13 AfQ group. Bottom right: Dwarf galaxy orbiting in a 10 13 Mq group. 



and b s fEa. lB2(l . However. lGill et all (|2004f) present peri- 
center distributions for subhalos at z = in cluster-sized 
dark matter halos with M v<gr = 1 - 3 x 10 14 Af Q . They 
find that the so distribution can be fit with a Gaussian 
with so — .35 and o~ So = .12. This distribution shows 
little variation with either M„ i9r or M VtSat . The max- 
imum ram pressure in a cluster with ICM parameters 
corresponding to a few 1O 14 M and c g = 7.5 increases 
by a factor of 5.4 when the pericenter decreases from 
s = 0.47 0.23. 

Some satellites have not yet passed through the peri- 
center of their orbits, and the maximum ram pressure 
they have experienced is determined by their current 
s or bit- Increasing s or &# from 0.35 to 1 results in a de- 
crease of the mass at which a satellite can be stripped by 
two orders of magnitude. Therefore, unless it is known 
that a group of satellites have all passed through their 
pericenters, the effect of varying s or bit cannot be treated 
like the scatter in the other parameters. Instead, when 
comparing satellites in different environments or with dif- 
ferent masses, s or bu must be specified. 

At a given s or uu the ram pre ssure is dependent on 
e/e v , but not l/l v . In simulations iVitvitska et all l)2002fl 
find that the speeds of satellites entering dark halos 
can be fit with a Maxwell-Boltzmann distribution with 
v rms ~ l-15u^, v 2 c = GM Vigr /r Vl g r . The total energy 
per unit mass of an object at the virial radius with 



this speed is e « e v , and is a weak function of c g . 
This distribution has a large scatter, a v = .6v c , corre- 
sponding to e/e v between 0.2 and 2.8. At s or bit = -35, 
when v sa t{s or bit = l)/«c is varied between 0.55 and 1.75, 
the M v S at /M y nr needed for stripp ing vary by a fac- 
tor of ~ 1.6. IVitv itska et al. (2002) also find that for 
M V}Sat < AM v<gr the angular momenta of satellites en- 
tering halos can be fit with a Maxwell-Boltzmann distri- 
bution with Irms/lv = -7, l v = v c r v . As expected, the 
distributions of e/e v and b s do not depend on M VtSat 
and M Vygr . However, unless the correlation between 
v S at(s rbit = ^)/ v c and l/l v is understood, these param- 
eters cannot be used to determine the distribution of so- 

4.3. Inclination 

The Gunn and Gott condition assumes that all galaxies 
orbit in a face-on wind. In reality, galaxies orbit at all 
inclinations. Simulations show that galaxies are stripping 
at all inclinations, but edge-on galaxies are stripped to 
l arger radii than f a ce-on galaxies. 

Marcolin i et al.1 l)2003ft point out that while galaxies 
that are hit face-on by the ICM wind experience the most 
ram pressure stripping, galaxies in an edge-on wind ex- 
perience the most viscous stripping. They derive an ex- 
pression for the ram pressure needed to viscously strip 
an edge-on galaxy to radius R and treat this pressure 
as a restoring pressure for the edge-on case. When they 
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compared the new condition to the Gunn and Gott con- 
dition, they find that any ram pressure that can com- 
pletely strip a face-on galaxy will viscously strip an edge- 
on galaxy and that any ram pressure that cannot strip 
a face-on galaxy will not strip an edge-on galaxy. The 
two models differ when a galaxy is partially stripped. In 
this case edge-on galaxies are still stripped, but to larger 
radii than face-on galaxies. They run simulations at a 
variety of inclination angles, i, and find that the final 
stripping radii match the predictions of the two condi- 
tions. When galaxies at intermediate angles are partially 
stripped, they are stripped to intermediate radii. 

When Marcolini et al. plot the two restoring pressures 
for three galaxy masses, 0.76, 7.4, and 77.2 xl0 9 M Q , 
their Figure 3, the fractional difference between the two 
restoring pressures is close to constant throughout the 
outer disk but increases with the mass of the galaxy. For 
the largest galaxy the restoring force in the edge-on case 
is « 4.5 times larger. 

4.4. Concentration of the Cluster 

Small dark matter halos ten d to be more concentr ated 
than large dark matter halos. iBullock et alJ l|2001f) use 
an analytical model of the evolution of dark matter halos 
to describe this dependency. For the standard ACDM 
model, near M v = M* w 1.5 x 10 13 ft -1 M Q they find 

* = 0) = 9 (t)" 3 (26) 

The c g used for the clust er and group models are deter- 
mined using this relation. IBullock et alJ l)2001fl also see a 
scatter that is as large as the evolution of the concentra- 
tion over the range O.OlAf* < M v < 100A/*. To explore 
the scatter introduced through c g , c g is varied between 5 
and 16. 

The cluster concentration affects both sq and the ram 
pressure at a given s or bit (eas. IB2I and IB6f) . When c g 
is increased from 5 to 16, for v1 at {s or bit = 1) = 1.15^ 
and l/l v = 0.7, the pericenter decreases from sq — 0.53 
to 0.5. This is smaller than the scatter in the pericen- 
ter. At a given s or bit, the orbital speed increases slightly 
as c g increases. This difference increases as s or bit de- 
creases. For Sorbit — 0.75, the change is not noticeable. 
For s or bit — 0.35, the M v ^ at /M Vi g r at which a satellite is 
stripped increases by a factor of 1.2 when c g is increased 
from 5 to 16. 

4.5. Concentration of the Satellite 

In the outer H I disk the restoring pressure of the 
dark matter halo makes a non-negligible contribution 
to the restoring pressure. More concentrated satellites 
have higher restoring pressures than other satellites of 
the same mass. Varying c s between 12 and 20 varies 
the M v sa t /M v gr at which stripping occurs by a factor 
of w 1.7. 

The concentration also alters how the hot galactic halo 
is stripped. In a more concentrated dark matter halo the 
gas is denser in the center and more diffuse elsewhere. 
Stripping can therefore occur down to a smaller radius, 
but a smaller fraction of the gas is lost. 

4.6. Length Ratios, Xi 



How disk scale lengths vary with M„ jSat is an open 
question. Theory suggests that when cold gas disks form 
in dark matter halos, their initial radii are pr oportional 
to th e virial radii of the halos they form in l)Mo et alJ 
1998(). There is a large scatter in this relationship that 
is due to a large scatter in the disk angular momenta. 
However, cold disks are composed of a combination of 
stars and gas, and how these components arise from the 
original cold disk is not understood. 

The fraction of the H I that is found in the outer disk, 
where the restoring force from the stellar disk is weakest, 
is set by the ratio Xd/X g (eqs. El and [TJJ| . iSwaters et alJ 
l|2002j) observed the H I disks of 73 local dwarfs. They 
found Xd/Xg = 6±2.5 for the dwarfs. For late-type spirals 
^d/X g = 5.5±1.6. These observations suggest that Xd/X g 
contributes to the scatter in M VtSa t/M Vtgr , but not to the 
mass dependence. Varying Xd/X g changes both the size 
and the shape of the restoring force. When A g is held 
at 5.7 and Xd/X g is decreased to 3.5, the restoring force 
decreases within R/R g = 0.45 and increases at larger R. 
The fractional change in P res t is less than 1.1 throughout 
the disk. Holding X g constant and increasing Xd/X g to 
8.5, results in decreasing the restoring force outside of 
R/Rg = 0.3 by a factor of less than 1.13. 

Varying X g while holding Xd/X g constant changes the 
density of both the gas and stellar disks. At a given 
R/Rg, the density of both disks cx A^ and the restoring 
pressure cx Xg . The scatter in X g and how X g evolves with 
mass are not known. However, because of the heavy de- 
pendence of the restoring pressure on it, the scatter in X g 
is likely to make an important contribution to the scatter 
in the M VtSa t/M Vtgr at which stripping occurs. Increas- 
ing A g from 0.75A S to 1.25A g decreases the M v ^ sa t/M v ^ gr 
needed for stripping by a factor of 20. 

The presence of molecular gas complicates the physics 
of stripping in the inner disk. Therefore, this model is 
only applied to disk radii beyond 1.5 stellar disk scale 
lengths. At these radii, R > lOrb, the bulge acts as a 
point mass, and changing A& (Eq. [SJ has no effect. In 
particular, the radius of the bulge changes with the mass 
of the bulge. Below, the mass of the bulge is allowed to 

1/3 

vary by 40%. Assuming that r& cx M b , this corresponds 
to allowing to vary by 12%. This is not noticeable in 
the outer H I disk. 

4.7. Mass Fractions, nii 

Both mb and md g feas. Inland II 1J) are correlated with 
the mass of the satellite. While the majority of large 
spiral galaxies have a significant fraction of their stars in 
a stellar bulge, low-mass late-type galaxies do not. Very 
low m ass spiral galaxies exist that do not contain any 
bulge ({Matthews fc Gallagherlll997]) and dwarf galaxies 
are observed to be either dE/dSph or dlrr. In addition, 
low-mass late- type galaxies have a higher f raction of their 
mass in the gas disk than large spirals ijSwaters et alJ 
2002). Decreasing mj reduces the restoring pressure in 
the outer disk. In contrast, increasing m,d g increases the 
restoring pressure. 

In the bottom right panels of Figures and [21 the 
satellite is bulgeless, the gas mass fraction, mdg, is dou- 
bled, and the stellar mass fraction, nrids, is reduced by 
20%. The original model corresponds to a large late- 
type spiral, while the bulgeless model is closer to a late- 
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type dwarf galaxy. Because the restoring force due to the 
bulge is small, the net effect is to increase the restoring 
pressure. Most galaxies will lie between the two plotted 
models. 

The n%i can be split between those that set the 
depth of the potential well, m<j s + mH 2 and mb, and 
rridg which sets the gas density. In late-type galax- 
ies, bulges can contribute up to 50% of the stellar 
light ijBinnev fc Mcrrificld 1998). In the reference model, 
25% of the stars are in the bulge. Varying mj, between 
15% and 35% of the satellite's stars, varies the restor- 
ing force in the outer disk by a factor of w 1.1. Neither 
the mass dependence nor the scatter in rrids + tuh 2 are 
known. Varying m,d s + itih 2 by 20%, varies the restor- 
ing force by a factor of ss 1.3. The effect of scatter in 
the gas density parameter, nid g , is straight forward. The 

M V}Sat /M v gr at which stripping occurs varies as m d ^ 2 . 
Varying mdg by a factor of 1.5, similar to the scatter 
in Mhi/Lr in Figure 9 of Swaters et. al., varies the 
M V:Sat /M Vy g r needed for stripping by a factor of 1.8. 

Variations in the m, values should be correlated be- 
cause the variation in the mass in baryons in a satellite 
is likely smaller than the variation in the mass of a given 
component. The correlation is such that it will reduce 
the overall scatter. 

5. DISCUSSION 

5.1. Comparison with Observations and Simulations 

In this section, the M V)Sa t/M Vtgr values at which the 
model predicts stripping sho uld occur are compa red to 
observations of H I disks from lSolanes et al] (J2001) . This 
sample is useful for making comparisons between the 
model and observations because the number of spirals 
is large and the observations span several clusters. The 
average temperature of the clusters is 3.5keV. This cor- 
responds to M V: g r = 3 x 10 14 Mq. The average gas 
profile for this temperature can be described using the 
middle mass cluster model. In Figure 4 from Solancs 
et al. the fraction of spirals with def > 0.3 and the 
average deficiency are plotted versus r/RA, where Ra 
is the Abel radius. In Figure 01 of the present paper 
M V}Sat is plotted versus the deficiency in a 3 x 10 14 M Q 
cluster at several s or bit- The deficiency is defined as 



def =_\og(M /M bs), where M a bs is the observed H I mass 
and M is the average H I mass for field spirals with the 
same morphology and optical diameter. For the model 
def = log(M g /(M g -M str )). 

The observed fraction of spirals with def > 0.3 reaches 
0.3 at r/R A « 1. and 0.45 at r/R A « 0.5. The model 
predicts that almost all galaxies within s or bu = 0.5 
have def > 0.3. However, only w 30% of the galax- 
ies viewed within r/r vgr = 0.5 actually reside within 
r / r v,gr = 0.5. This estimate assumes that the galaxy 
distribution follows an NFW profile and that all galaxies 
within r/r v , gr — 2 are included in the projected clus- 
ter. At SorUt = 1, the model predicts that an aver- 
age galaxy traveling face-on to the wind is stripped to 
def > 0.3 for M VtSat < lO n - 5 M . This estimate can be 
roughly adjusted for inclination by multiplying by 10~' 5 , 
M VtSat < 10 11 M Q . Approximately 45% of the galaxies 
seen within s or bit — 1 actually reside within s or bu = 1- 
Some galaxies beyond the s or bit under consideration will 
be stripped to def > 0.3. Therefore, the observed frac- 
tion of spirals with def > 0.3 should be greater than 
that predicted by multiplying the volume fraction with 
the fraction of stripped galaxies within s or bit- This con- 
dition is met for s or bit = 1 if at least a third of the ob- 
served spirals have M v>sa t > 10 11 Af Q , corresponding to 
Md ~ 8 x 10 9 AT Q . For s or ut = 0.5, the volume fraction 
is lower than the observed fraction and this condition is 
easily met. The difference between the observed value 
and product of the volume fraction and stripped frac- 
tion should decrease as s or bit increases. A IO^Mq spiral 
traveling face-on to the wind is stripped to def w 1 at 
Sorbit = 0.8. A 10 12 M Q spiral traveling face-on to the 
wind is stripped to def rs 1 at s n rb it = 0-4. Spira l s with 
def k, 1 are not observed in the iSolanes et all (2001) 
sample outside of IRa- 

The model does a good job of matching the observa- 
tions of H I disks from lSolanes et all l|200lD . If the model 
over-predicted the degree of stripping, then it would over- 
predict the fraction of spirals with def > 0.3 and predict 
deficiencies at large radii greater than those observed. If 
the model under-predicted the degree of stripping, then 
it might not predict that spirals in the center of the clus- 
ter can have def ~ 1 and would require that a larger 
fraction of the observed spirals reside in the cluster cen- 
ter. 

Bot h iMarcolini et all l)2003|) and iRoediger fc Henslerl 
( 2005) run simulations of disk galaxies in winds typi- 
cal of the outskirts of a group, equivalent to s or bit — 1 
and 0.7, and the center of a group, s or bit = 0.2 and 
0.25. In Marcolini et al. the two lowest mass simu- 
lated dwarfs are completely stripped in the high ram 
pressure. In the other cases, the dwarfs are partially 
stripped to a range of radii. When the gas surface den- 
sity in the current model is adjusted to the exponential 
gas density in the dwarf models, the current model and 
the simulations match well for both s or bit ■ The only dis- 
parity is that for s or bit = 0.2 the high-mass model is 
stripped to R s tr/Rg — 0.1 in the simulations and to 
R s tr/Rg — 0.45 in the model. This is probably due 
to differing dark matter halo profiles. The high-mass, 
M > 10 u Mq, galaxy model of Roedinger and Hensler is 
stripped to R s tr/Rg = 0.8 in the lower ram pressure and 
Rstr/Rg = 0.6 in the higher ram pressure. After the gas 
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density is adjusted, the stripping radii of their high-mass 
model in the two winds are matched well by the model. 

5.2. Model Predictions 

In this model there exist three stripping regimes, strip- 
ping of the hot galactic halo, the flat H I disk, and the 
inner gas disk. The focus is on stripping of the flat H 
I disk. However, the model also includes a hot galactic 
halo. This hot gas is easily stripped in all environments. 
A satellite galaxy that is orbiting at s or ut — 1- and 
has a hot galactic halo consisting of 0.5% of its mass, is 
stripped of at least 40% of its hot halo. At s or ut = 0.35, 
a typical pericenter, all satellite galaxies are stripped of 
practically their entire halo in both groups and clusters. 
The stripping of the inner disk has not been studied. 
However, this gas should be difficult to strip and any 
spiral that has not been stripped of its entire outer H 
I disk should retain this gas. The rest of this section 
focuses on the outer disk. 

The model does a reasonable job of matching H I ob- 
servations in large cluster spirals. Therefore, it can be 
used to predict the extent of stripping at different M VySat 
and M Vy g r . The model's basic prediction, that the strip- 
ping radius is determined by the ratio M VySat /M Vygr is 
modified by the mass dependence of the model's param- 
eters. Therefore, the usefulness of the model's predic- 
tions depends in part on how well this dependence is 
understood. Fortunately, the mass dependence of many 
of the parameters is based on observations. The model 
parameters that most alter M s tr{M VySat /M Vy g r ) in low- 
mass systems are the ICM profile parameters and the 
distribution of mass in the galaxy. The scaled ram pres- 
sure, Pram I M^gr i that a satellite experiences decreases 
with M Vy g r because the ICM density throughout most 
of the cluster decreases. The scaled restoring pressure, 

P res t/M^s at , is larger in low-mass disks because a higher 
fraction of the disk mass is found in the H I disk. The 
changes in both the scaled ram and restoring pressures 
act to decrease the M VySat /M Vy g r at which stripping oc- 
curs in low-mass systems. For example, the Local Group 
dwarfs have M VySat / M VyQr smaller than the bright spi- 
rals in Virgo but should experience similar amounts of 
stripping. 

The model makes three basic predictions. Dwarf galax- 
ies are stripped of their entire outer H I disks in clus- 
ters and lose varying fractions of their outer disk in 
groups. Massive spirals can also be stripped of signifi- 
cant fractions of their H I disk in groups if they travel 
to small s or bit- In a 10 13 M Q group, a 10 10 M Q dwarf is 
stripped of 50% of its H I disk at s orhit = 0.7, a IO n M 
satellite is similarly stripped at either s or bit — 0.5 or 
Sorbit = 0.35, depending on the satellite model used, 
and a large 5 x 10 11 M Q satellite is similarly stripped 
at s rbit — 0.3. The corresponding orbits from stripping 
80% of the H I disk are s or Mt = 0.55, 0.35, 0.25, and 0.2. 
It should be kept in mind that these numbers are for an 
average satellite traveling face-on to the ICM wind. 

The degree of stripping that a particular satellite of 
mass M v sat experiences orbiting in a cluster of mass 
M v . gr can vary greatly. The parameters that are respon- 
sible for most of the variation in M str are the ICM pa- 
rameters, a, /?, r Vy g r /r c , the galaxy's inclination, i, the 
pericenter, sq, and the extent of the disk, X g . The varia- 



tion in the stellar and gas surface densities is largely con- 
tained in X g . The ICM parameters are set by the choice 
of cluster, but clusters contain galaxies with a variety of 
orbits, morphologies, and inclinations. Within the same 
cluster, a galaxy of mass M VySat = M with a low surface 
density on a radial orbit traveling face-on to the wind can 
be as severely stripped as a galaxy with M VySat ~ M /2b 
that has a high surface density and is traveling edge-on. 
At the same time, the ram pressure can vary significantly 
between clusters with the same mass because of varia- 
tions in the ICM profile. Two identical satellite galaxies 
on identical orbits in different clusters of the same mass 
can experience ram pressures that vary by a factor of 
w 10. 

An effective mass, M e ff, can be defined for each in- 
dividual galaxy. This is the mass of a satellite that 
the model predicts is stripped to the same R s tr as the 
given galaxy. The effect of the scatter in the param- 
eters can then be discussed in terms of the scatter in 
M e ff at a given M VySat . The scatter in the M VySat /M Vygr 
at which a galaxy is stripped of M str , or in M e ff, can 
best be expressed as a multiplicative factor, log(M e ff) — 
log(M VySat /M Vy g r ) ± oiog. The source of scatter in M e ff 
can be grouped into three types, that due to the ICM 
parameters, a, /?, r Vygr /r c , the orbital parameters, e/e v , 
i, c g , and the restoring pressure parameters, c s , mi, 
Xi. The overall scatter in the high- and low- mass mod- 
els, the result of varying all of the model parameters, 
are <Ji og ,h ~ 0.75 and 07 O9 ,a w 0.85. These can be 
broken down into aiogjCMM = 0.7, criog,iCM,i = 0.8, 
o'log.orbit = 0.5, and o-io gy rest = 0.65. The discussion in 
the paragraph above is based on these estimates. 

The uses of this model are twofold. First, it makes gen- 
eral predictions about the degree of stripping that occurs 
to galaxies of different masses in a range of environments. 
Second, it links the degree of stripping, M str or R str , to 
the ratio M VySat /M Vy g r . This may be useful for verifying 
whether stripping is occurring within a galaxy catalog. 
The size of the scatter in M e ff determines the size of 
the galaxy catalog that is needed to see a dependence of 
M str on M VySat /M Vygr . The scatter in the (3 profile pa- 
rameters is a major contributor to the scatter in M e ff. 
Therefore, such a catalog must either contain many clus- 
ters, or the ICM profile of the clusters must be known. 
As an estimate of the catalog size necessary, if 200 satel- 
lite galaxies with approximately the same M VySat orbiting 
in 50 groups or clusters with the same M Vygr are com- 
bined, then log(M e ff) w log(M VySa t) ± -25. If the range 
in M VySat and M Vy9r of such a catalog is large compared 
to this scatter, then the effects of stripping should be 
seen in the relationship between M str and M VySat / 'M Vygr . 

6. CONCLUSIONS 

The model developed here relates the degree of ram 
pressure stripping a satellite galaxy experiences to the 
galaxy and cluster masses and can be used to quickly 
determine the extent to which a galaxy is likely to be 
stripped. All galaxies lose most of their hot galac- 
tic halo. In clusters galaxies at moderate mass ra- 
tios, M VySat /M Vy gr w 10 -2 , are moderately stripped, 
Rstr/Rd ~ 0-6, at intermediate distances, 0.5 < s or bit < 
1, from the cluster center and severely stripped closer to 
the cluster center. Satellites with lower M VySat /M Vygr are 
severely stripped even at intermediate distances. Strip- 
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ping also occurs in groups. However, the same degree of 
stripping occurs for lower M VtS at/M v>gT in groups than 
in clusters. Dwarf galaxies are moderately to severely 
stripped at intermediate distances in groups, and large 
spiral galaxies can be moderately stripped if they travel 
to small s or bit- 

The model is simple and motivated by observations. 
Dark matter profiles are well matched to NFW profiles 
outside a possible core, observations of X-ray gas in clus- 
ters can be fit using (3 profiles, and the average gravita- 
tional potential of most galaxies should be well matched 
by the model potential. The model has a large num- 
ber of parameters, but the values for many can be taken 
from observations. However, the model assumes that 
clusters are static. From an evolutionary standpoint, 
it is only valid after the group or cluster has acquired 
an ICM. On the other hand, in dynamic clusters bulk 
ICM motions may cause more stripping to occur than is 
predict ed by the model. Th is is observed in the Virgo 
cluster ijKennev et alJl2004b|) . 

Ram pressure stripping is not the only way to 
remove gas from galaxies. In groups and clus- 
ters ti dal stripping of bo t h gas and stars can occ ur 
(e.g., ijBureau et all 120041: IFatterson fc Thuanl \199^ ). 
For dwarf galaxies it has been proposed that supernova 
winds a ssociated with bursts of star formation may ex- 
pel g as ijFerrara fc Tolstovl feOOO: Sili ch fc Tenorio-Taglel 
|2001|) . and the efficie ncy of this mechanism may depend 
on the environment l)Murakami fc Babull 11993) . How- 
ever, ram pressure stripping is capable of removing gas 
from galaxies across a large range of galaxy masses and 
environments. In particular ram pressure stripping can 
act when supernova-driven ejection is likely to be ineffi- 
cient and can remove gas from galaxies that are either 
not tidally interacting or that are experiencing only weak 
tidal interactions. 

The tidal radius, r t , at the physical orbital radius, 
r or bit, as estimated using r t = (M ViSat /3M enc )r or bit, can 
be compared to the stripping radius, R s t r , as deter- 
mined using the model. Here M enc is the cluster mass 
enclosed within r or ut- In scaled coordinates, r t /Rd = 
^gSorbit[3m(s or bit)]~ 1/3 , where m(s or ut) = M enc /M Vigr . 
The H I fraction lost to tidal stripping depends only on 
the scaled size of the gas disk and s or bit and not on ei- 
ther mass. Because the gas disk resides in the center 
of the satellite's dark halo, tidal stripping rarely effects 
the H I disk and practically never competes with the ef- 
fect RPS. For example, in order for the H I disk to be 
tidally stripped to r t /R g — .8, a satellite must travel to 
Sorbit ~ -1. By this s or bit the entire outer H I disk has 
been ram pressure stripped for almost all satellites. 

The effect of varying the model's parameters was stud- 
ied both to identify the parameters that most effect the 
extent to which an individual galaxy is stripped and to 
determine the range of satellite and cluster masses that 
result in the same R s tr- The extent to which an in- 
dividual galaxy is stripped depends most strongly on 
the galaxy's s or Mt or sq, inclination, i, and disk scale 
length, X g , and on the cluster's ICM profile. Galaxies 
that reach smaller s or bit, are face-on to the wind, have 
denser disks, or orbit through a denser ICM are more 
effectively stripped. Galaxies in the same cluster with 
M V}Sat that differ by as much as a factor of 25 can be 



stripped to the same R s t r - In different clusters, galaxies 
on identical orbits with identical morphologies, identical 
mi, Ai, c s , can be stripped to the same R s tr with M VtSat 
values that differ by as much as a factor of 30. 

Ram pressure stripping of satellite galaxies' outer gas 
disks and hot galactic halos is occurring frequently in 
both groups and clusters. In general, removing gas from 
a galaxy reduces star formation. The gas in the outer 
disk and hot galactic halo is not currently involved in 
star formation, but may feed future star formation if it is 
not stripped. Therefore, it would be interesting to study 
in more detail how the removal of this gas affects galaxy 
evolution and to determine if, to what extent, and how 
quickly star formation declines after a galaxy is stripped. 
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APPENDIX 
A. NFW PROFILE 

When simulations are run that allow cold, non-interacting particles to evolve gravitationally from a flat power 
spectrum at high redshift, the dark matter halos that emerge at low z are well fit by a universal density profile over 
many orders of magnitude in mass. This profile was originally parameterized bv [Navarro et all l|1997j) . Both the 
satellite and group/cluster dark matter halos are described using an NFW profile, the relevant properties of which are 
summarized here. 

The NFW profile is 

P{r) _ $char 

p° (r/r s ){l + r/r s y 1 ' 

where p° c is the average density of the universe at z = 0, r s is an inner scale radius, and 5 c h a r is a characteristic 
over-density. It can be rewritten in terms of a concentration, c, an overdensity, v p , and a scaled radius, s = r/r v 



p{s) _ v p c 2 g(c) 



p° c 3s(l + cs) 2 
c = — 



(A2) 
(A3) 



v p c 3 g{c) 

Ochar — ^ V ) 

5(c) = , . 1 J7T—, (A5) 
ln(l + c) - c/(l + c) 
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Hester 



The overdensity, virial mass, M v , and virial radius, r v , are related by 

M v = -m*v p pl 



Note that the virial radius scales as M v ' . The mass contained within s is 

cs 



r l/3 

M(s) = g(c)M v 



ln( 1 + cs) - 



and the gravitational potential is 



*(*) 



1 + cs 
GM v g{c) ln(l + cs) 



The above can be found in ijCole fc LacevHl996t ILokas fc Mamonll200H) and ijNavarro et al.lll997|) . 

B. ORBITS IN AN NFW POTENTIAL 

Bound orbits in a central force field travel between two radial extremes. These are found by solving. 

1 _ 2 [e - $(r)] 
^2 - p 

where e and I are the energy and angular momentum per unit mass l)Binnev fc Tremaindll988l) . 
The scaled radial extremes, sq, for an NFW profile are given by 



1 _ -2 



where e v and b s are defined as 



e v ln(l + cs ) 

1 g(c) 

e s 



GM V 



\e\rl 



The parameter b would be the impact paramete r for an unbound orbit. 
Conserving angular momentum and using Eq. IB2I the orbital speed at the pericenter, is 

2 _ b 2 s e GM V 
'7> — 



J — 2 



1 /3 

It is proportional to M v and depends on the concentration through sq(c) 
At any point in the orbit, the orbital speed, v sa t, is 



2GM V 



5(c) 



ln(l + cs) 



C. HOT GALACTIC HALO 
The virial temperature, T v , of a gravitational potential is defined by 



Sk B T v 



\W\ 



(A6) 

(A7) 
(A8) 



(Bl) 

(B2) 

(B3) 
(B4) 

(B5) 
(B6) 



(CI) 



where W is the work done by gravity in forming the halo, Ur is the Boltzmann constant, m p is the mass of the proton, 
and /i is the mean molecular weig ht ijBinnev fc Tremaindll988|) . 
Using eqs. IA2I IA6I and lA7l the virial temperature of an NFW profile is 



\W\ _ GM V 2 2 



M v 



A dimensionless temperature, t, is defined. 



cV(c) 



1 ln(l + cs) - cs/(l + cs) 



s(l + cs) 2 



ids 



GAL, 



v 2 



t(c) 



GAL 



t(c) = 



c 2 g(c)h(c) 



(C2) 
(C3) 

(C4) 
(C5) 
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TABLE CI 

Parameter Values and Scatter 



pjireiiiiGtGr 


ref in text 


reference model 


■A A 1 

mid mass model 


1 A 1 

low mass model 


scatter 


AP a 


AM a 


a 


Eq.@ 




6.5 


5.5 


20 


x6 


6 


15 


B h 


Eq.@ 




0.67 






±.08 


2 


3 












0.57 


_ 








_ 








_ 




0.45 


±.08 


7.4 


20 


1 b 


Eq.@ 




12 


25 




7-20 


6 


16 














170 


75 - 270 


5.4 


13 




sec. 




0° 






0° - 90° 


4.5 


10 


Vsat(s or bit — l)/ u c 


Eq. 




3 


1.15 


_ 


_ 


±.6 


1.4 


1.6 


C 9 


Eq. 






5.5 


7.5 


10 


5-16 


1.1 


1.2 


c s 


Eq. 






12 




20 


12-20 


1.4 


1.7 




Eq. 






5.7 






±1.4 


8 


20 




Eq. 






6 






±2.5 


1.2 


1.4 




Eq. 






0.04 












x b 


Eq. 






180 




na 










Eq. 






0.008 




0.016 


xl.5 


1.5 


1.8 


m ds + m H2 


Eq. 






0.05 




0.04 


±.01 


1.3 


1.1 


m b 


Eq. 






0.017 







0.012-0.024 


1.1 


1.2 


m 3 g 


Eq. 






0.005 






xl.5 


1.5 


1.8 


so 


Eq. 


B2l 


0.35 






±.12 


5,1 
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a AP and AM columns list the scatter in P rarn or P res t and M v ^ sa t/M v _ gr . The scatter is given as the fractional change when 
the parameter of interest is varied over its entire range. Therefore, log(AM) £S 2(T; 03 (M U)Sat /M t , iflr ). 
^The first row shows affects of scatter for the large mass model and the third row for the small mass model. Columns 7 and 8 



show the fractional increase in P r , 



at s or bit 



For a gas in hydrostatic equilibrium in potential 0, ignoring self-gravity 

PM=PoCxpfe[0(r)-0(O)] 



ln(l + cs) 



Combining eqs. IC3l and lC6l the gas density profile is 



p g (s) = p exp 



ch(c) 



Letting the gas mass within r v equal a fraction m sg of the virial mass results in a central density of 



Po 



Pb v p m sg 
3/ 2 (c) 



The integral, /2(c), is 



12(c) = J s 2 dsexp 
The scaled density profile j is defined by 

1 

j(s,c) 



ch{c) 
3 



h(c) 6XP Vc/i(c) 



ln(l + cs) 



In (cs) 



cs 



P( s ) = ( Vj ^r ) m sgj(s,c) 



(C6) 

(C7) 

(C8) 
(C9) 

(CIO) 
(Cll) 



